How to generate and measure anomalous weakly non-ergodic 
Brownian motion in simple systems 

A. Fulinski 

M. Smoluchowski Institute of Physics, Jagiellonian University, 
Reymonta 4, PL-30-059 Krakow, Poland 

Abstract 

It is shown that in systems with time-dependent and/or spatially nonuniform temperature 
T(t,x), (i) most of the transport processes is weakly non-ergodic, and (ii) the diffusion (Brow- 
nian motion, BM) is anomalous. A few examples of simple arrangements, easy for experimental 
realization, are discussed in detail. Proposed measurements will enable also the observation of 
transitions from ergodic to weakly non-ergodic and from normal to anomalous diffusion. New 
effects are predicted: (i) zero-mean oscillations of T(t) accelerate BM (pumping effect), (ii) the 
combination of temporal and spatial variations of temperature may lead to superballistic BM, (iii) 
linear gradients of T{x) result in an exponential acceleration of BM. One can expect similar effects 
in inflationary systems with time-dependent metrics. 
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In the recent literature 



3| the process X(t) is called weakly non-ergodic when: 



lim (S*X{t,t f )?(X\t)), (1) 



tf—tOO 



(thorouhout this paper X(t — 0) = 0) with 

8 2 X(t, tf) = — !— - ^[X( s + t)-X(,)] 2 (2) 
tf-t J 

relating two positions of the process (walker) separated by a time lag t (for more details cf. 

3). 

Weak ergodicity breaking (WEB) was discussed mainly in the context of anomalous diffu- 
sion (stochastic theory) 2] , though recently it was shown that WEB is typical for much wider 
class of cumulative-growth processes, both stochastic and deterministic [3] . Experimentally, 
WEB was found in composite disordered systems: various glasses, cytoplasm of living cells, 
and like [4]. In this note we want to show how both anomalous diffusion (Brownian motion) 
and WEB can be generated and measured also in much simpler systems, easier to deal with 
and better suitable for various experimental manipulations. 

Consider Brownian motion (BM) in a medium with variable temperature. The thermal 
(white Gaussian) noise driving the Brownian particle is characterized by temperature- 
dependent intensity: 

(£(*)> = 0, mm = o* e (T)5(t - s) . (3) 

Consider first simple setups which can be easily realized experimentally: (i) thin (quasi- 
2d) flat layer (film) of a fluid which is cooling or heating (e. g. by emitting or absorbing 
radiation), (ii) long thin rod or fluid-filled cylinder immersed in a heat bath [6]. When 
the characteristic time for equilibration of momenta is much shorter than the characteristic 
hydrodynamic (diffusional) time, it is safe to assume that the thermal local equilibrium is 
established with well-defined local temperature (in a coarse-grained sense). Note that when 
the cooling/heating is spatially uniform, there will appear no temperature gradients through 
the system, and we shall get the system in a quasi-adiabatically shifting thermal equilibrium. 

In the considered system cr| will be changing in time through T = T(t). For an ideal 
medium cr| ~ T (fluctuation-dissipation theorem), and we may write cr|(T) = T/T <7q(To), 
where To > is some reference temperature, e. g., T(T = 0); however, in real systems such 
dependence can be more complicated, depending on details of interparticle interactions. In 
a spatially non-uniform system the noise intensity will also depend on the (coarse-grained) 



position r, = a^(r,t) (time-dependent diffusion coefficient was used recently for the 
description of random walks in porous materials [3]). 

The Brownian motion driven by the noise with = cr^(t) is weakly non-ergodic. Namely, 
consider kinetic (Langevin) equation jlj]: 

± s X p (t)=t^ (t), > -1/2, (4) 

where £ is the thermal noise at constant reference temperature T , with dispersion <Tq(T ). 
The formal solution of Eq. (J3J) is 

s X (t) = [ u%{u)du. (5) 



Jo 

This scaled Brownian motion (SBM) is equivalent to the Brownian motion with time- 
dependent diffusion coefficient |5j: define the scaled noise rj(t) = t^^ (t). Then 

( V (t)) = , ( V (t)v(s)) = o*{t)6(t - s) , a v (t) = ^+V , (6) 

and 
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j tS x,{t) = n(t) . (7) 



It was shown in [3[ that SBM is weakly non-ergodic. This proof can be easily extended for 
any Brownian motion driven by thermal noise with time-dependent intensity c^(t) = 
4>(t)a (in fact this follows directly from the formulas (2.2)-(2.3) of j^j]; for the sake of 
completeness, we present here the specialized formulas). Namely, for 



X(t)= / ds<j>(s)£ { s), (8) 
Jo 

we get from Eqs. p, (J2U and ([S]): 

(X 2 (t)) = 2 fdsD{s), (9) 



o 



2D(t) = a&it) = al(t,T(t)) } (10) 



tf-t ps + t 



{5 2 X(t : t f )) = / ds / duD(u) 

if — t J0 Js 

1- f * ds[{X\s + t))-{X\s))}, (11) 
f f ~t Jo 
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which is weakly non-ergodic except for D = const, i. e., for the Wiener process (with 
constant temperature), and maybe for some special combinations of the spatial and temporal 
dependence of D(r,t) (i. e., T(r,t) - cf. Eq. (123|) below with 7 = 2/3 for an example). 

When the heating function D(t) is proportional to $ (SBM case), this setup enables 
the realization of both WEB, and anomalous diffusion. Other forms of <f>(t) will lead to 
non-diffusional (i. e., with dispersion not scaling algebraically with time) behavior of the 
(weakly nonergodic) Brownian motions. 

Interesting is the case of periodic changes of temperature. Let <f>(t) = 1 + Asin(cut), 
A 2 < 1. Then 

A A 
(X 2 {t)) = 2D A t + AD - [1 - cos(wt) - - sin(2wt)l , (12) 

ui 4 

where Da = D (l + A 2 /2), 2D = ctq(T ), which tends asymptotically to normal diffusion 
with higher effective diffusion coefficient (enhanced diffusion or accelerated Brownian mo- 
tion). The same effect will appear for anomalous Brownian motions, and for other forms of 
driving <f>(t), either periodic or random. Note that when the asymptotic temporal average 
of the changes of temperature is zero, as in Eq. fll2|) . we have a kind of pumping of energy 
into the Brownian particle, analogous to e. g. ionic nanopumping in asymmetric nanopores 



More challenging for the analysis (though perhaps experimentally simpler) is the case 
when T = T(r,t). Analytical treatment is possible for the separable T, when T(r,t) = 
T (j) 2 : (t)ip 2 : (r) , and only when the appropriate stochastic differential (Langevin) equation: 

dX(t)/dt = 0(t)V(r)£ o (t) (13) 



is treated according to the Stratonovich interpretation [9] (other situations need to be treated 
either by numerical simulations or by numerical solutions of the appropriate Smoluchowski- 



Fokker-Planck equations 10, llll])- In that case, the Eq.( fl3l) can be written in the form: 



dY{t)/dt = <p(t)i{t) , dY{t) = dX/ip(r) , (14) 

which means that the process {Y(t)} is equivalent to the scaled Brownian motion, i. e., 
is weakly nonergodic 3|. To obtain properties of the original process {X(t)}, inversion of 
the transform (fT4l) is necessary. Note that for stationary systems, = const., the process 
{Y(t)} becomes the standard Wiener process {W^(£)}, which is ergodic. This, however, does 
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not imply that all original stationary processes {X(t)} are also ergodic. Namely, the formal 
solution of Eq.f lT3|) reads: 



X(t)= / ds*(;(X(s))cl>(s)Z(s). (15) 
Jo 

Making use of the Stratonovich interpretation of stochastic integrals [9|, we get 



(X 2 (t)}=2D / ds^iXis))}^), (16) 
Jo 



which is an integral equation for the dispersion of X. 

Analogous calculations of the time average, Eq.([2]), give: 



{px(t,t f )) = r ^([x(,+t)-x( S )] 2 ) 

l f - 1 Jo 



2D , f s+t 



ds I du(tlj 2 (X(u))}(j) 2 (u) . (17) 

tf ~t Jo Js 

Thus, even for the stationary case = 1, the process {X(t)} is ergodic only when the second 
integral in Eq. ffTT|) ^ f{s), which is fulfilled for if) = const., but besides maybe only for some 
very special choices of the temperature gradients. This proves that in most realistic cases the 
processes {X(t)} in the media with spatial temperature gradients are weakly non-ergodic. 

Further analysis depends on the forms of the functions 4>(t) and - especially - ^(r). Note 
that, when D depends on r and t through its dependence on temperature: D = D(T(r,t)), 
the position- and time-dependent temperature profile in a heated/cooled system has to be 
calculated from the Fourier-Kirchhoff equation (FKE): 

dT d 2 T . 

_ = A — + S( rii ), (18) 



where A is the heat conduction coefficient and S - heat source [12j. The solution of this 
equation for an arbitrary shape of the source S can be found only by numerical computations, 
especially when we want to heat/cool the system in a very narrow region in the immediate 
neighbourhood of r = 0. On the other hand, the FKE can be used for an inverse problem: 
one may calculate from Eq. (1181) the space- and time-dependent heat flow S(r,t) necessary 
for the attaining the prescribed temperature profile. In principle, any required T(r,t) can 
be used. However, while the calculation of the appropriate heat flow is a simple task, its 
experimental realization might not be so easy. 

As an example, consider a simple quasi-one- dimensional system: heat-conducting long 
narrow rod (or a fluid contained in a narrow cylinder) of length 2L (x G [— L, +L], L — > oo), 



with initial condition T(X = 0,t = 0) = To = const, and X(t = 0) = 0. To retain 
the correspondence with anomalous diffusion, assume that T = T o 2 (t)-?/> 2 (X) = T (t + 
t) 2/3 (|X| + a) 27 , r > 0, a > 0. Then a? 1 1$ = 1, and T is positive and bounded. 
From Eq.flHD: 

Y(t) = sgn(X(t))\ (19) 
[ln(|X(t)|/a+l), 7 = 1. 

The inversion to X(t) gives: 

(sgn(Y(t)/n)l\Y(t)/^-a},, 
= I (20) 
{sgn(X{t))a[e Y ^ s ^(^ - l] , 7 = 1 

where /i = 1/(1 — 7 ), and sgn(f) denotes the sign of / (sgn(0) = 0). Note that the initial 
condition X(t = 0) = implies Y(t = 0) = 0. 

n 

The results for scaled Brownian motion mentioned above [3] give 

(Y 2 (t)} = 2Dp(t + r) a ° , « = 2/3 + 1 , (21) 

and for tf — > oo, t <C tf 



(5*Y(t, t f )) -> 2D p {t + r) {t f + r)— 1 . (22) 



where Z)g = D /a . Thus [13], for 7 7^ 1 

(X 2 (t)> = 2D M (t + to) a , 

a = /i« = (2/3 + l)/(l-7) 



2D M = {2D p /n 2 Y r((2/i + l)/2)^V, (23) 



and for 7 = 1 



(X 2 ) = exp{4 J D^(t + r) ao }- 

- 2aexp{D /3 (t + r) a °} + a 2 . (24) 

Note that 7 = 1 implies parabolic increase of temperature T(X) with distance from the 
center. 



On the other hand, direct calculations of 5 2 X(t,tf)} cannot be done analytically [ijj]. 
However, to prove that the original process {X(t)} is weakly non-ergodic, it is sufficient to 
prove that the ensemble and time averages of any function of this process are not equal. 
Therefore the result ( |2Tl) -(j22l) for the averages of Y(t) = f(X(T)) proves the weak ergod- 
icity breaking of diffusional processes in systems with spatially nonuniform and/or time- 
dependent temperature T(t,x). Note however that for 7 = —2/3 we get a — 1, i. e., the 
process {X(t)} becomes normal diffusion. 

The results discussed in this paper show that the anomalous nonergodic transport is 
not restricted to composite systems, but can be easily generated in almost any system. 
Therefore, by (uniform) heating of the system with changes of temperature proportional to 
t 13 one can measure directly the weakly-nonergodic anomalous diffusion (both subdiffusion for 
— 1/2 < /3 < 0, slowing-down the heating, and super diffusion for < (3 < 1/2, accelerating 
the heating). 

A corollary to this implies that the Mandelbrot- Van Ness (cf. js|) fractional Brownian 
motion Bfjit) which is ergodic, will become weakly non-ergodic in systems with varying 



temperature. Viz., the MBM of [3| can be written as Bh driven by f](t), Eq.(E) 



B H (t) 



1 



(t - u) H -^ 2 ri(t]T(t))dt- 



T(H + 1/2) 


(-uf- 1/2 r](t;T(t))dt 



00 

t 



T(H + 1/2) 


(-uf- 1/2 <P(u)dW(u 



{t~u) H ^ 2 <j){u)dW{u)- 



(25) 

where dW{t) = £(t,T )dt, and W(t) is the Wiener process at temperature T . 

Consider now a few specific situations, which lead to new, so far unknown effects, which 
will appear in systems with variable temperature. 

(i) First such effect - the enhancement of diffusion due to temperature oscillations was 
discussed above (Eq.f fT2"]) )). The spatial analog of this situation, i. e., the stationary system 
with periodic (or aperiodic) variations of temperature, can be realized experimentally in the 
so-called bistable ballast resistor, i. e., thin conducting wire with nonlinear dependence of 
resistivity on temperature. When such systems are heated by electric current, and cooled 
by surface radiation, there appears the transition to the so-called dissipative structure with 



interlaced hot and cold fragments 15] . Unfortunately, the analytic treatment of the Brow- 



nian motion in such situations (i. e., the inversion of the series {Y(t)} into {X(t)}, and 
subsequent averaging over the noise) is practically impossible. 

(ii) Another effect is the superballistic Brownian motion resulting from the cooperation 
of the spatial and temporal changes of temperature. From the discussed above case of 
anomalous diffusion, when <p(t) = (t + r) 13 and ip{X) = (\X\ +a) 7 we find (cf. Eq. (I23j) ) that 
a > 2 when + 7 > 1/2. 

(iii) The results described by Eq.f l23|) show that the Brownian motion accelerates ex- 
ponentially in stationary systems with parabolic increase of temperature from the center 
towards the ends of a long (formally, infinitely long) cylinder. This means also that the 
logarithm of the dispersion of X(t) behaves asymptotically as the normal diffusion. Linear 
T(X) = T (|X| + a) results in the parabolic acceleration of the Brownian motion (ballistic 
motion) . 

As we have mentioned, the analytic treatment of the Brownian motion in systems with 
spatial temperature gradients is in most cases very difficult (or just impossible). However, 
although the analytic calculations of dispersions (both ensemble and temporal) are too 
difficult, it is possible to find the averages of some functions of X(t), e. g. the inversions of 
the process Y(t), or some functions of these inversions. 

All these assertions can be checked experimentally. This can be done by standard mea- 
surements of the time-dependence of the dispersion (X 2 (t)) and - when possible - by the 
measurements of the single trajectories of Brownian particles. The latter would give di- 
rect experimental verification of non-ergodic character of SBM, and of other processes of 
this type. Note that the same can be done with other time characteristics of the speed 
of heating/ cooling, including natural exponential cooling by radiation. Note also that it is 
important (and obvious) in such measurements that the time-scales of heating/ cooling and 
of diffusion should be of the same order of magnitude. 

It was shown here that anomalous diffusion (and related to it weak ergodicity breaking) 
can be detected in simple systems with nonstationary or/and nonuniform temperature. It 
should be noted that (i) similar effect can be found also in systems driven by external 



forces 



16|, (ii) time-dependent diffusion coefficient was used recently for the description of 



random walks in porous materials |7f, (iii) one can expect similar effects in structures with 
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time-dependent metrics (inflationary systems). 
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